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Abstract 

Governing equations are derived for the kinetics of physical aging in polymeric 
glasses. An amorphous polymer is treated as an ensemble of cooperatively rearranged 
regions (CRR). Any CRR is thought of as a string of elementary clusters (EC). Frag- 
mentation of the string may occur at random time at any border between ECs. Two 
string can aggregate at random time to produce a new string. The processes of aggre- 
gation and fragmentation are treated as thermally activated, and the rate of fragmenta- 
tion is assumed to grow with temperature more rapidly than that for coalescence. This 
implies that only elementary clusters are stable at the glass transition temperature T g , 
whereas below T g , CRRs containing several ECs remain stable as well. A nonlinear 
differential equation is developed for the distribution of CRRs with various numbers of 
ECs. Adjustable parameters of the model are found by fitting experimental data for 
polycarbonate, poly (methyl methacrylate) , polystyrene and poly (vinyl acetate). For 
all materials, fair agreement is established between observations and results of nu- 
merical simulation. For PVAc, the relaxation spectrum found by matching data in 
a calorimetric test is successfully employed to predict experimental data in a shear 
relaxation test. 

1 Introduction 

This paper is concerned with the kinetics of structural relaxation (physical aging) in amor- 
phous glassy polymers. Slow dynamics in out-of-equilibrium disordered media [supercooled 
liquids, structural (including polymeric) glasses, disordered ferromagnets and antiferromag- 
nets, orientational glasses, vortices in superconductors, dipolar glasses, liquid crystalline 
colloids, spin glasses, etc.] has attracted substantial attention in the past decade, see, e.g., 
surveys 0, |3|, f| and the references therein. Unlike most previous studies which focused 
on equilibrium thermodynamics of glass-like systems ||, the present work deals with the 
evolution of physical properties of amorphous polymers quenched from some temperature 
T above the glass transition temperature T g to a temperature T in the sub-T g region and 
isothermally annealed at the temperature T. Changes in the mechanical response, specific 
volume, configurational entropy and other features of glassy polymers caused by thermal 
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jumps have been widely studied (both experimentally and theoretically) in the past four 
decades. We would like to mention here seminal works by Kovacs || [7| and Struik ||. For 
a detailed review, the reader is referred to |J. 

Despite a number of publications concerned with the evolution of internal structure of 
disordered media after thermal jumps, see @ for a survey, it is difficult to mention a model 
(phenomenological or molecular) which correctly predicts experimental data (for example, 
memory and chaos phenomena in cyclic thermal tests |10|, |TT| ) and establishes an adequate 
correspondence between observations in various (mechanical, dielectric, calorimetric, dilato- 



metric, etc.) experiments, see, e.g., [12|, [L3|, [TJ], [15], [16], [17) and the references therein. 

An important class of constitutive models for structural relaxation in amorphous mate- 
rials employs the concept of coarsening in disordered systems First models belonging 
to this class (the so-called droplet models) have been proposed a decade ago by Koper and 
Hilhorst [T8[] and Fisher and Huse JB| to predict the response of spin glasses. According to 
their approach, an arbitrary macro-region of a disordered medium is treated as an ensemble 
of droplets (domains with non-regular boundaries possessing non-zero fractal dimensions) 
where spins are aligned with one of the ground states. These regions alter with time because 



of slow coarsening, and the characteristic size of a region grows as a power of time |L8 



or as a logarithm of time [19 . For a detailed review of phenomenological models for the 



phase-ordering kinetics, we refer to . 



Another approach to the description of domain growth has been developed by Ben-Nairn 
and Krapivsky |H], ^4]]. The present study employs basic ideas of their works to 

predict enthalpy recovery in amorphous polymers and to study the effect of waiting time t w 
on the material response in mechanical (relaxation) tests. We would like to stress two issues 
that distinguish our constitutive equations from those derived by Ben-Nairn and Krapivsky: 

1. In the analysis of coalescence, we presume that for any domain, the rate of aggregation 
with another domain is proportional to the probability density (not to the concentra- 
tion) of regions with a certain volume (length). This makes solutions of constitutive 
equations independent of the total number of elementary clusters at the initial instant. 

2. We combine the aggregation-fragmentation theory with the concept of internal time. 
This allows fair approximation of experimental data for the relaxing enthalpy to be 
established (the conventional approach fails to adequately describe enthalpy recovery 
during 4 to 5 decades of time). 

The concept of internal time is widely used to describe thermo-mechanical response of poly- 
mers (for a review, the reader is referred to [p5fl). In this study, we employ a model with an 



entropy-driven material clock. This approach was proposed by Struik || and was successfully 
used to predict shear-thickening of polymer solutions in |26| . 



With reference to the theory of cooperative relaxation J27j, an amorphous polymer is 
thought of as an ensemble of cooperatively rearranging regions (CRR) that relax at random 
times as they are thermally agitated. A CRR is thought of as a globule consisting of scores 
of strands of long chains |23 . The characteristic length of a CRR in the vicinity of the glass 
transition temperature amounts to several nanometers pSfl . 

We introduce a minimal number of strands which permits rearrangement, u, and suppose 
that the number of strands in any CRR reads nu, where n is an integer and v is the number 
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of strands in an elementary cluster (EC). The concept of elementary domains was first 



introduced by Adam and Gibbs |27| for amorphous polymers and by Kadanoff |30] for spin 
glasses. For a discussion of the Kadanoff conjecture (near the critical point, large blocks of 



spins behave like individual spins) see [[31], |32|]. Description of the aging process in glassy 



polymers in terms of models with finite numbers of states was suggested by Chow |33] and 



Robertson [[34 



Unlike previous studies where the domain growth was interpreted as random motion of 
their boundaries annihilated at the contact points, see, e.g., f24j| , we treat structural recovery 



as a result of two (thermally activated) processes with different dependences of their rates 
on temperature: fragmentation of large regions and aggregation of small ones. An indirect 
confirmation of this picture is provided by experimental data in calorimetric tests at various 
temperatures Tk in the sub-T g region. Observations reveal that if T2 > 7\, the relaxing 
enthalpy AH in a test with the temperature T\ exceeds that in a test with the temperature 
T2 at small times, and the inverse is true at large times. 

The exposition is organized as follows. In Section 2, constitutive equations are derived 
for enthalpy recovery in amorphous glassy polymers. These equations are verified in Sec- 
tion 3 by comparison with experimental data for polycarbonate, polystyrene, poly(methyl 
methacrylate) and poly(vinyl acetate). The purpose of fitting observations is two-fold: (i) 
to confirm that the model can correctly predict measurements, and (ii) to analyze the effect 
of temperature T on material parameters. In Section 4, the model is subject to a more 
elaborate evaluation: after determining adjustable parameters by matching data in calori- 
metric tests, we use the probability density of CRRs with various energies to fit the material 
response in a static mechanical test. To make this possible, a model is derived for stress 
relaxation in disordered media based on the traps concept (for detail, we refer to 
Some concluding remarks are formulated in Section 5. 



2 Constitutive equations for domain growth 

The following hypotheses are introduced to develop a master equation for the concentration 
of CRRs: 

1. Structural relaxation in a disordered medium is governed by two processes at the 
micro-level: fragmentation of CRRs and their aggregation. 

2. Denote by P(t,n) the number of CRRs (per unit mass) at time t consisting of n + 1 
elementary clusters. The quantity P(t, n) is the mass concentration of CRRs containing 
n borders between ECs subject to fragmentation (to simplify the analysis, any CRR 
is thought of as a linear "string" of ECs). The function P(t, n) obeys the conservation 
law 

00 

X>+l)P(i,7l)=So, (1) 
n=0 

where S is the number of ECs per unit mass. 

3. The fragmentation process is characterized by its rate 7 that equals the number of 
fragmentation acts per boundary between ECs per unit time. We assume 7 to be a 
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function of the current temperature T only. The rate of changes in the quantity P(t, n) 
induced by fragmentation is given by 



Vf(t,n) = j[-nP(t,n) + 2 £ P{t,m)]. (2) 

m=ra+l 

The first term on the right-hand side of Eq. (2) determines the number of CRRs 
(containing n boundaries) destroyed by fragmentation, whereas the other term is the 
rate of creation of CRRs with n boundaries caused by fragmentation of relaxing regions 
containing larger number of elementary clusters. The physical meaning of formula (2) 
is discussed in detail in 



4. Coalescence of relaxing regions occurs when CRRs (with n and m borders between 
ECs) merge and create a new CRR (with n + m + 1 borders). The rate of aggregation 
of an individual CRR with n boundaries between ECs with CRRs with m borders 
is proportional to the density of CRRs with m boundaries (the ratio ip(t,m) of the 
number of CRRs with m boundaries to the total number of CRRs per unit mass) and 
to the number of random meetings, T njTn , for two CRRs with n and m boundaries. The 
rate of changes in the quantity P(t, n) induced by coalescence reads 

oo 

V c (t,n) = -P(t,n)[ J2 ^n, m <p(t,m) + 2Y n>n ip(t,n) 

m=0,mj^n 

n-1 

+ r mi „_ m _iP(t, m)(p(t, n - m - 1) (3) 



m=0 



with 



P(t,n) 

The first two terms on the right-hand side of Eq. (3) determine the number of CRRs 
destroyed by aggregation with other CRRs (the coefficient 2 in the second term means 
that two relaxing regions with n borders disappear when they meet one another). The 
last term in Eq. (3) characterizes the rate of creation of CRRs with n boundaries by 
aggregation of two relaxing regions with smaller numbers of ECs. 

5. The rate T n , m by which regions with n and m boundaries aggregate decays exponen- 
tially with the growth of the indices n and m, 

r n , m = Z-7exp -X(n + m) , (5) 

where L and A are material parameters. Assumption (5) is conventional for the treat- 
ment of (short-range) interactions between CRRs, but it differs from standard relations 



based on the theory of anomalous diffusion, see, e.g., |39| and the references therein, 
which lead to the power law dependence of T m ^ n on n and m. Our numerical analysis 
shows that replacement of Eq. (5) by a power law does not affect significantly the 
accuracy of fitting. 
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6. We employ the concept of material (reduced) time and suppose that 7 is the rate 
with respect to some internal time r. The rate of fragmentation with respect to the 
"universal" time t reads 

7o = 07, (6) 

where a is a shift factor. 

7. For polymers with an entropy-driven material clock ||, we set 

In a = —KqS, (7) 
where k is some material parameter and s is the configurational entropy per EC. 

8. The configurational entropy s is defined by the Boltzmann's formula 

00 

s(t) = -k B 53 n) In n), (8) 

n=0 

where k B is Boltzmann's constant. It is worth noting some difference between Eq. (8) 
and conventional relations, see, e.g., PD |. In Eq. (8) the configurational entropy s 
is defined per elementary cluster, whereas traditional formulas define configurational 
entropy per CRR. 

It follows from Eqs. (2) to (7) that the current concentration of CRRs (per unit mass) 
is governed by the differential equation 



dP 



-nP(t,n)+2 P^™; p ,. 

m=n+l l^m=0 r l r i m 



X 



m=n+l 



P(t,n) P{t,m)exp(-\(n + m)^ + 2P 2 (t, ra) exp(-2Ar 

m=0,mj^n 

n—1 -1 n 

-exp(-A(n-l)) ^2P(t,m)P{t,n-m-l) \, (9) 

where k = IcbKo- 

The study is confined to one-step "quench-and-wait" tests, when a polymer equilibrated 
at some temperature T > T g is quenched to a temperature T < T g and is annealed at the 
temperature T, 

m = { T *\H do) 

Assuming that above the glass transition temperature T g , only elementary clusters may exist, 
we postulate 

P(0, n) = S 5„, , (11) 

where 5 n ^ m is the Kronecker delta. The solution P(t,n) of Eqs. (8), (9) and (11) is indepen- 
dent of the total number of ECs, S , and it is determined by 4 adjustable parameters: L, 
To, k and A. 
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Let h(t) be the configurational enthalpy per EC. Assuming h to be expressed in terms of 
the configurational entropy s by the conventional relationship, 

— = T, 

ds 

and integrating this equality for the thermal program (10), we find that 

h(t) = Ts(t) {t > 0). (12) 

Multiplying Eq. (12) by the concentration of ECs, Ho, we calculate the configurational 
enthalpy per unit mass. Assuming the relaxing enthalpy AH (measured in calorimetric 
tests) to coincide with the configurational enthalpy, we arrive at the formula 

oo 

AH(t) = -A ^ ip(t, n) In (p(t, n) (13) 

n=0 

with 

A = k B TE . (14) 



3 Comparison with observations in calorimetric tests 

We begin with experimental data for polycarbonate. For a detailed description of samples 
and the experimental procedure, see []4Tj] . Adjustable parameters are found by matching 
observations using the steepest-descent procedure. Figure 1 demonstrates that the model 
correctly predicts experimental data at two temperatures in the sub-T g region. Setting 
T = T g = 420 K in Eq. (14), using mass density p = 1.196 g/cm 3 []4"2"|| , and the value A = 0.8 
J/g found by fitting observations, we obtain 

H = 1.15 • 10 26 m~ 3 . (15) 



This value is in accord with H = 3.6 • 10 26 m~ 3 for polytetrafluoroethylene [4"3"|, but it 
is less than concentrations of holes measured for a series of polycarbonates using positron 
lifetime spectroscopy |44| (for comparison, we assume that any EC may be associated with 



a micro- hole). To explain this discrepancy, we recall that (i) Bohlen et al. [44j presumed a 



Gaussian distribution of holes with a peak far above the zero volume, whereas our calculations 
demonstrate that the probability density of CRRs substantially differs from the Gaussian 
ansatz (see Figure 2), and (ii) the value (15) is for temperatures in the vicinity of T g , while 
the PALS measurements were carried out at room temperature. 

Suppose that the effect of temperature on the rate of fragmentation 70 is described by 
the Arrhenius formula 

70 CO =7*exp(-— ), (16) 

where AE is some activation energy and R is the universal gas constant. Applying Eq. (16) 
to two temperatures, T\ and T 2 , in the sub-T g region, we find that 

A B =^L>„2^. (17) 

T 2 - Ti 70 
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Taking values of 7o(T]t) determined by fitting experimental data (Figure 1), we obtain AE = 
174.1 kcal/mol, which is in fair agreement with AE = 173.6 kcal/mol found by shift of 
relaxation curves measured at various temperatures in the vicinity of the glass transition 
temperature [p5] . 

According to Figure 1, an increase in temperature T leads to a sharp decrease in the 
relative rate of aggregation L and to an increase in the parameter A. This conclusion is 
in fair agreement with our hypothesis that above the glass transition temperature T g only 
elementary clusters exist (because the rate of merging vanishes, the growth of CRRs caused 
by their coalescence is prohibited). It is worth noting that an increase in A with temperature 
observed in experiments makes questionable models of aggregation that are based on the 
diffusion mechanism, see ||39|| , provided that the coefficient A in Eq. (5) may be associated 
(at least, to some extent) with the decrease in the diffusion coefficient induced by the growth 
of mass of domains to be aggregated. 

The equilibrium distributions depicted in Figure 2 are characterized by their mean value 
Mi and variance M 2 , 

00 00 
M x = 5^71^(00, n), M 2 = 5^(n-Afi)V(°°,n)- (18) 

n=0 n=0 

Numerical simulation demonstrates that the quantities Mx and M 2 substantially increase 
with a decrease in T. This is in accord with the conventional approach to the description 
of structural relaxation [fj, |5j which presumes a substantial growth in the roughness of the 
energy landscape with a decrease in temperature (which, in turn, results in a monotonic 
increase in the variance of the distribution of CRRs with AT = T g — T) . 

We proceed with matching experimental data for polystyrene using the above algorithm. 
For a description of the experimental procedure, we refer to f45fl . Figure 3 reveals fair 



agreement between observations at two temperatures in the sub-T g region and results of 
numerical analysis. Equation (14) with T = T g = 373.0 K, p = 1.04 g/cm 3 and A = 1.1 J/g 
results in So = 2.55 • 10 26 m -3 , which is in acceptable agreement with PALS measurements. 
It follows from Eq. (17) and Figure 3 that the activation energy AE = 136.4 kcal/mol, 
which is rather close to the value AE = 170.5 kcal/mol found by shift of creep curves in the 



sub-T g region [46]. The equilibrium distributions of CRRs containing various numbers of 
elementary clusters are plotted in Figure 4 which demonstrate that an increase in AT leads 
to an increase in Mi and M 2 . This is in qualitative agreement with conventional models 
for aging that predict the growth of roughness of the energy landscape with a decrease in 
temperature. 

Experimental data for polycarbonate and polystyrene (only at two temperatures in the 
sub-Tg region) are not sufficient to analyze the effect of temperature on the kinetics of struc- 
tural relaxation. To study changes in adjustable parameters caused by the annealing tem- 
perature T, we match observations for poly (methyl methacrylate). For a detailed description 
of samples and the experimental procedure, we refer to p7| . We begin approximation with 
measurements at the lowest temperature T = 375.0 K and determine parameters L, 70, k, A 
and A using a version of the steepest-descent algorithm. Afterwards, we fix A (an analog of 
the number of elementary clusters per unit mass) and fit data at other temperatures using 4 
constants: L, 70, k and A. The assumption that the parameter k is temperature-dependent 
differs the treatment of data for poly(methyl methacrylate) from that for polycarbonate and 
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polystyrene (for which k was assumed to be independent of temperature). This is done 
to ensure better quality of matching experimental data. Figures 5 to 10 demonstrate fair 
agreement between observations and results of numerical simulation. 

The rate of fragmentation 70 is plotted versus the increment of temperature AT in 
Figure 11. This figure shows that experimental data are correctly approximated by the 
"linear" function 

log7o = a - aiAT (19) 

with adjustable parameters a^. Equations (16) and (19) imply that in the vicinity of the 
glass transition temperature T g , the activation energy AE is given by 

AE = ai RTl In 10. 

It follows from this formula and Figure 11 that AE = 70.81 kcal/mol, which is in excellent 
agreement with data for activation energy obtained by shift of creep and relaxation curves at 
various temperatures. For example, using WLF parameters for poly(methyl methacrylate) 
48] , we arrive at the value AE = 83.2 kcal/mol, which is close to our result. This is 



in contrast with the Cowie-Ferguson model [47| which results in the value AE = 164.8 
kcal/mol for the same set of experimental data. 

The parameter L (which characterizes the ratio of the rate of aggregation to the rate of 
fragmentation) is depicted in Figure 11 versus the increment of temperature AT. Experi- 
mental data are fairly well approximated by the linear function 

L = b + b x AT (20) 

with adjustable parameters bk- Equation (20) implies that L monotonically increases with a 
decrease in temperature T, which means that the rate of aggregation of elementary clusters 
grows with a departure from the glass transition temperature. 

The dimensionless parameters A and k are plotted versus the increment of temperature 
AT in Figure 12. Experimental data are approximated by the functions 

A = c - ci AT, log k = C - C]AT, (21) 

where and are adjustable parameters. Figure 12 demonstrates that the quantities A 
and k increase with temperature. Far below the glass transition point (about 30 K below 
T g ) the parameter A vanishes, which results in an essential simplification of Eq. (9) (whose 
coefficients become independent of n and m, which allows a method of generating functions to 
be applied to develop analytical solutions [f49"|). Explicit solutions to the governing equations 



are, however, beyond the scope of the present study. 

It is worth noting that phenomenological models for physical aging in polymers presume 
that the effects of temperature and structure may be separated, see, e.g., f7j. Results of 
numerical simulation reveal that this hypothesis is valid when the increment of temperature 
AT is not too small (say, it exceeds 15 K) which implies that k becomes weakly dependent 
on temperature. 

Conventional constitutive equations for structural relaxation in amorphous polymers em- 
ploy the random energy model |50|| , according to which the distribution of relaxing regions 



with various energies is Gaussian with temperature-dependent mean value and standard 
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deviation. Fitting experimental data obtained in mechanical tests for a number of glassy 
polymers show that the standard deviation of energy of CRRs, S, linearly decreases with 
temperature and vanishes at some critical temperature T cr slightly above T g |37] . Because the 
present study is based on the assumption that only elementary clusters are stable at the glass 
transition temperature (which implies that T cr = T g ), our purpose now is to demonstrate 
that 

S = \[M~2 

linearly increases with AT in the sub-T g region. Figure 13 shows that with an acceptable 
level of accuracy the graphs Mi (AT) and E(AT) are linear. Approximation of experimental 
data by the functions 

Mi = mi AT, S = m 2 AT (22) 

yields m 2 = 0.15 K _1 . This result is in good agreement with 777,2 = 0.09 K _1 for polycarbonate 
and 7772 = 0.53 K _1 for poly(vinyl acetate) found by fitting observations in static and dynamic 
mechanical tests fl37| . 

Hitherto, only average characteristics (first moments and entropy) of the distribution of 
CRRs have been compared with observations. To reveal that the distribution of relaxing 
regions tp(t,ri) itself is in agreement with measurements, we analyze experimental data for 
poly(vinyl acetate) in calorimetric and mechanical (relaxation) tests. For a description of 



the experimental procedure, see |L6 



We begin with matching observations in a calorimetric test at T = 303 K. Figure 14 
demonstrates fair correspondence between experimental data and results of numerical sim- 
ulation. Setting T = T g = 315 K in Eq. (14) and using the values p = 1.182 g/cm 3 []5T| and 
A = 1.6 J/g (Figure 14), we arrive at S = 3.24 • 10 26 m~ 3 , which is in good accord with the 



concentration of holes measured by PALS for polytetrafluoroethylene ||43|| . All adjustable 



parameters for PVAc (except for the rate of fragmentation 70) are similar to those for PC, PS 
and PMMA. The quantity 70 is smaller than that for other polymers. This may be explained 
by the differences in the glass transition temperatures of materials under consideration (glass 
transition points for PC, PS and PMMA substantially exceed T g for PVAc) and the fact that 
fragmentation and aggregation of CRRs are thermally activated. 

To fit data in mechanical tests, stress-strain relations should be derived for amorphous 
glassy polymers under short-term mechanical loading (compared to the waiting time t w ). In 
this work, we briefly sketch the development of constitutive equations confining ourselves to 
uniaxial relaxation tests. Greater detail of derivations can be found in 107 . 



4 The mechanical response in relaxation tests 



In accord with the hoping concept, see, e.g., |2"%|] , a relaxing region is treated as a point 
located at the bottom level of its potential well on the energy landscape. Rearrangement of 
CRRs is modeled as hops of relaxing regions to higher energy levels. Hops occur at random 
times, and they are driven by thermal fluctuations. Below the glass transition temperature, 
energy barriers between potential wells are assumed to be so high that thermally agitated 
CRRs cannot leave their traps. Adopting the transition-state theory |p2[ , we suppose that 
rearrangements occur when CRRs reach some liquid-like (reference) energy level. The posi- 
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tion of this level is not fixed, but it slowly ascends in time approaching some limiting value 
as t w — > oo. 

Any potential well is described by its depth w > with respect to the position of the 
reference energy level at the initial instant (immediately after the quench). Following our 
treatment of CRRs as aggregates composed of integer numbers of elementary clusters, we 
assume that the set of available energies w is countable, 

w = wx,w 2 , ■ ■ ■ ,w n , . . . , 

and the value w n is proportional to the volume (length) of an appropriate CRR (string), 

w n = an, (23) 

where a > is an adjustable parameter. 

Let q{uj)duj be the probability for a CRR to reach (in a hop) the energy level that exceeds 
its bottom level by an energy to 1 located in the interval [u, u + du}. Referring to the extreme 
value statistics 0, we set 

q(uj) = Aexp(-Auj), 

where A is a material constant. The probability to reach the reference state in a hop for a 
CRR trapped in a potential well with the depth w reads 

f-OO 



Q(t w ,w) = / q(uj)duj = exp — A[w + £l(t 



/u>+fi(t w ) 

Here fi(t w ) is the increment of the reference energy level after the waiting time t w with respect 
to that at the initial instant, and we take into account that the duration of conventional 
mechanical tests is negligible compared to t w . 

The rate of hops in a potential well, T, is defined as the number of hops (of an arbitrary 
intensity) per unit time. Assuming T to be independent of w (according to the theory 
of thermally activated processes, T is a function of the current temperature T only), and 
multiplying the rate of hops T by the probability of reaching the liquid-like state in a hop 
Q(t w ,w), we arrive at the Eyring formula for the rate of rearrangement |53" |, 



R(t w ,w) = T (t w )exp(-Aw), T = T expL4fi(t w ) ) . (24) 



The rate of rearrangement, R(t w ,w), can be thought of as the ratio of the number of rear- 
ranging regions (per unit time) to the number of CRRs to be rearranged, 

R(t w , w) = -— -J -^(t, t w , w). (25) 

For relaxation tests, the quantity E(t,t w ,w) is the number of CRRs (per unit mass) located 
in cages with energy w that have not been rearranged until time t (which is measured from 
the beginning of the test). Integration of Eq. (25) implies that 

E(t,t w ,w) = H(0, t w , w) exp —R(t w ,w)t . (26) 
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The initial condition for Eq. (26) is given by 



E{0,U,w n ) = N(t w )p(t e ,n), (27) 

where 

oo 

iV(t w ) = ^P(t w ,n) 

n=0 

is the number of CRRs (per unit mass) after annealing for the time t w . 

The natural configuration of a CRR after rearrangement is assumed to coincide with the 
deformed configuration of the viscoelastic medium at the instant of rearrangement. This 
implies that rearranged CRRs are stress-free in a relaxation test, and only non-rearranged 
regions have non-zero mechanical energies. The strain energy density of an amorphous 
polymer (per unit mass) under uniaxial loading, U, equals the sum of mechanical energies 
of individual regions. A CRR is thought of as a linear elastic medium with the potential 
energy of deformation 

1 2 

where \i is the rigidity and e is the macro-strain (for definiteness, shear deformation is 
considered). Combining this expression with Eqs. (23), (24), (26) and (27), we find that 

oo 

U(t,t w ) = -fie 2 J2 E(t,t w ,w n ) 

* n=l 

■y OO 

= -fiN(t w )e 2 n) exp -T (t w ) exp(-a n)t , (28) 

1 n=l 

where cxq = a A. The stress a is expressed in terms of the mechanical energy U by the 
formula 

dU 

° = P ^- 

This equality together with Eq. (28) implies that 

oo 

G(t, t w ) = G (t w ) J2 V(*w: n ) ex P [- r o(^w) exp(-a n)t\ , (29) 

n=l 

where G(t, t w ) = cr(t)/e is the current shear modulus, Go(^w) = fiN(t w ) is the initial modulus. 

To verify the model, we calculate the distribution of CRRs, <p(t w ,n), using Eqs. (8), 
(9) and (11) with adjustable parameters found by fitting observations in the calorimetric 
test, and match experimental data in the relaxation test by Eq. (29). Given a constant Go, 
the material parameters r and «o are calculated using the steepest-descent procedure. The 
initial shear modulus Go is determined by the least-square algorithm. Figure 15 demonstrates 
fair agreement between observations and results of numerical simulation. 

It is worth noting that our constitutive equations employ two independent time-scales for 
the description of structural recovery in polymeric glasses. One scale is entirely determined 
by the distribution function (p(t w ,n). The other time-scale characterizes the ascent of the 
liquid-like energy level with respect to the energy landscape, and it is portrayed by the 
function Q(t w ). This approach is in agreement with that recently proposed by Tanaka |]5^}| , 
where two order parameters were used to predict slow dynamics in supercooled liquids. 
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5 Conclusions 



A model is derived for structural relaxation in amorphous glassy polymers after thermal 
jumps. A polymeric glass is treated as an ensemble of cooperatively rearranging regions 
whose concentration changes with time because of their fragmentation and aggregation. A 
CRR is modeled as a string (linear chain) of elementary clusters. Fragmentation of the 
string may occur at random time at any border between elementary clusters with equal 
probability. Aggregation of relaxing regions occurs at random time as well. The rate of 
coalescence for two CRRs decreases exponentially with the growth of their sizes. With a 
decrease in temperature T, the rates of fragmentation and aggregation decrease, but the 
rate of fragmentation reduces more rapidly. This implies that only elementary clusters exist 
at the glass transition temperature T g , whereas in the sub-T g region, CRRs consisting of 
several ECs may be stable as well. 

To verify constitutive equations, we fit experimental data for relaxing enthalpy for poly- 
carbonate, polystyrene, poly (methyl methacrylate) and poly (vinyl acetate). Fair agreement 
is demonstrated between observations in calorimetric tests and results of numerical analysis. 
Material parameters found by fitting measurements are in good accord with those determined 
experimentally in other tests. 

To establish correspondence between observations in calorimetric and mechanical tests, 
we find material parameters for PVAc by fitting relaxing enthalpy, determine the relaxation 
spectrum and use this spectrum to match data in mechanical (static) test. An acceptable 
agreement between experimental data for the shear modulus and numerical predictions con- 
firms our believe that the model may be employed for the analysis of physical aging in tests, 
where several experimental methods are applied simultaneously. 
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List of figures 



Figure 1: The relaxation enthalpy AH J/g versus time t h for polycarbonate annealed at the 



temperature T K. Circles: experimental data [[|TJ. Solid lines: numerical simulation with 
A = 0.8 J/g and « = 1.0. Curve 1: T = 408.0, L = 80.0, A = 0.3, 70 = 10 3 h _1 ; curve 2: 
T = 413.0, L = 5.0, A = 1.2, 7o = 1.2 ■ 10 4 br 1 

Figure 2: The equilibrium distribution of CRRs <p = <p (00, ri) for polycarbonate. Unfilled 
circles: T = 408 K, M x = 4.9328, M 2 = 9.8743. Filled circles: T = 413 K, M a = 0.9157, 
M 2 = 0.6758 

Figure 3: The relaxation enthalpy AH J/g versus time t h for polystyrene annealed at the 
temperature T K. Circles: experimental data [ 45] . Solid lines: numerical simulation with 



A = 1.1 J/g and « = 1.0. Curve 1: T = 363.0, L = 80.0, A = 0.3, 70 = 10 3 h" 1 ; curve 2: 
T = 366.0, L = 35.0, A = 0.5, 70 = 4.4 • 10 3 h" 1 

Figure 4: The equilibrium distribution of CRRs <p = <f(oo,n) for polystyrene. Unfilled 
circles: T = 363 K, Mi = 4.9328, M 2 = 9.8743. Filled circles: T = 366 K, M x = 2.7906, 
M2 = 3.6336 The relaxation enthalpy AH J/g versus time t h for poly(methyl methacrylate) 



at T = 375.0 K (unfilled circles: experimental data [57]; solid line: numerical simulation) 
and the final distribution of CRRs (p = <p(t w ,n) with t w = 10 2 ' 5 h 

Figure 5: The relaxation enthalpy AH J/g versus time t h for poly(methyl methacrylate) at 
T = 375.0 K (unfilled circles: experimental data [E7|; solid line: numerical simulation) and 
the final distribution of CRRs <p = <p(t w ,n) with t w = 10 2 5 h 

Figure 6:The relaxation enthalpy AH J/g versus time t h for poly(methyl methacrylate) at 
T = 377.5 K (unfilled circles: experimental data [57|; solid line: numerical simulation) and 



the final distribution of CRRs = <p(t w ,n) with t w = 10 2 ' 5 h 

Figure 7: The relaxation enthalpy AH J/g versus time t h for poly(methyl methacrylate) at 



T = 380.0 K (unfilled circles: experimental data [|7|; solid line: numerical simulation) and 



the final distribution of CRRs <p = <p(t w , n) with t w = 10 2 ' 5 h 

Figure 8: The relaxation enthalpy AH J/g versus time t h for poly(methyl methacrylate) at 



T = 382.5 K (unfilled circles: experimental data [ 57[ ; solid line: numerical simulation) and 



the final distribution of CRRs <p = <p(t w ,n) with t w = 10 2 ' 5 h 

Figure 9: The relaxation enthalpy AH J/g versus time t h for poly(methyl methacrylate) at 



T = 385.0 K (unfilled circles: experimental data ||7|]; solid line: numerical simulation) and 



the final distribution of CRRs if = (p(t w , n) with t w = 10 2 ' 5 h 

Figure 10: The relaxation enthalpy AH J/g versus time t h for poly(methyl methacrylate) 
at T = 387.5 K (unfilled circles: experimental data [57]; solid line: numerical simulation) 



and the final distribution of CRRs if = (p(t w ,n) with t w = 10 2 ' 5 h 

Figure 11: The rate of fragmentation 7o h _1 (unfilled circles) and the dimensionless param- 
eter L (filled circles) versus the increment of temperature AT K for poly(methyl methacry- 



late). Symbols: treatment of observations [47]. Solid lines: approximation of the experi 



mental data by Eqs. (19) and (20) with a = 4.3733, a x = 0.0993 and b = 3.8507 and 
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h = 1.2686 



Figure 12: The dimensionless parameters A (unfilled circles) and k (filled circles) versus 
the increment of temperature AT K for poly(methyl methacrylate). Symbols: treatment 
47 1 - Solid lines: approximation of the experimental data by Eq. (21) with 
0.7173, Cx = 0.0310 



of observations 
c = 0.7611, Ci 



0.0274 and C n 



Figure 13: The average number of borders between ECs in a CRR, M 1; and its standard 

i 

deviation, £ = M 2 2 , versus the increment of temperature AT K for poly(methyl methacry- 



Solid lines: approximation of the experimental 
= 0.1533. Curve 1: Mi; curve 2: E 



late). Circles: treatment of observations [f47 |. 
data by Eq. (22) with mi = 0.2037 and m 2 : 

Figure 14: The relaxation enthalpy AH J/g versus time t h for poly(vinyl acetate) annealed 
at T = 303 K. Unfilled circles: experimental data |16| . Solid line: numerical simulation with 
A = 1.6 J/g, k = 2.9, L = 20.0, A = 0.45, 7 = 4.74 ■ 10 2 h" 1 . Filled circles: the distribution 
of CRRs tp = (p(t w , n) with t w = 10 2 - 5 h, M x = 2.4209, M 2 = 3.2730 

Figure 15: The shear modulus G GPa versus time t s for poly(vinyl acetate) at T = 303 K 
and t w = 16.5 h. Circles: experimental data |li|. Solid line: prediction of the model with 
a = 2.24, T = 0.23 s" 1 and G = 1.0176 GPa 
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